Bose-Einstein condensation of atoms in a uniform potential 
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We have observed Bose-Einstein condensation of an atomic gas in the (quasi-)uniform three-dimensional 
potential of an optical box trap. Condensation is seen in the bimodal momentum distribution and the anisotropic 
time-of-flight expansion of the condensate. The critical temperature agrees with the theoretical prediction for 
a uniform Bose gas. The momentum distribution of our non-condensed quantum-degenerate gas is also clearly 
distinct from the conventional case of a harmonically trapped sample and close to the expected distribution in 
a uniform system. We confirm the coherence of our condensate in a matter-wave interference experiment. Our 
experiments open many new possibilities for fundamental studies of many -body physics. 

PACS numbers: 03.75.Hh, 67.85.-d 



Ultracold Bose and Fermi atomic gases are widely and suc- 
cessfully used as testbeds of fundamental many-body physics 
m . Experimental tools such as Feshbach interaction reso- 
nances f2l, optical lattices |3 |, and synthetic gauge fields 
offer great flexibility for studies of outstanding problems 
arising in many areas, most commonly in condensed matter 
physics. However, an important difference between "conven- 
tional" many-body systems and ultracold gases is that the for- 
mer are usually spatially uniform while the latter are tradition- 
ally produced in harmonic traps with no translational symme- 
tries. 

Various methods have been developed to overcome this 
problem and extract uniform- system properties from a har- 
monically trapped sample liSj-JJI, relying on the local den- 
sity approximation tSllTTIl or selective probing of a small cen- 
tral portion of the cloud ifTTilTSl . Sometimes harmonic trap- 
ping can even be advantageous, allowing simultaneous map- 
ping of uniform- system properties at different (local) particle 
densities. On the other hand, in many important situations lo- 
cal approaches are inherently limiting, for example for stud- 
ies of critical behaviour with diverging correlation lengths. 
The possibility to directly study a spatially uniform quantum- 
degenerate gas has thus remained an important experimental 
challenge. So far, atomic Bose-Einstein condensates (BECs) 
have been loaded into linear 1 14 | or circular 1 15] traps which 
are uniform along only one direction while still being har- 
monic along the other two. 

Here we demonstrate Bose-Einstein condensation of an 
atomic gas in a three-dimensional (quasi-)uniform potential. 
We load an optical box trap depicted in Fig. [TJa) with ^^Rb 
atoms pre-cooled in a harmonic trap and achieve condensa- 
tion by evaporative cooling in the box potential. Below a 
critical temperature, Tc ~ 90 nK, condensation is seen in 
the emergence of a bimodal momentum distribution and the 
anisotropic time-of-flight (TOE) expansion of the BEC. We 
characterise the flatness of our box potential and show that 
both the momentum distribution of the non-condensed com- 
ponent and the thermodynamics of condensation are close to 
the theoretical expectations for a uniform system, while being 
clearly distinct from the conventional case of a harmonically 
trapped gas. We also demonstrate a simple experimental con- 



figuration for trapped-atom interferometry and use it to con- 
firm the coherence of our quasi-uniform BEC. 

Our setup for producing ^^Rb condensates in a harmonic 
potential is described elsewhere ifTFIl : we create BECs in 
the |F, mi?) = |2,2) hyperfine ground state using a hybrid 
magnetic-optical trap ifTVl . The dark optical trap QSl [191 
which is central to this work is formed by three 532 nm laser 
beams - a "tube" beam propagating along the x axis and two 




FIG. 1: Preparing a quasi-uniform Bose gas. (a) The optical-box trap 
is formed by one hollow "tube" beam and two "sheet" beams creating 
a repulsive potential for the atoms. The atomic cloud, depicted in 
red, is confined to the dark cylindrical region. Gravitational force 
is cancelled by a magnetic field gradient B'. (b) The three trapping 
beams are created by reflecting a single Gaussian beam off a phase- 
imprinting spatial light modulator, (c) The atoms are loaded into 
the box trap after pre-cooling in a harmonic trap, (d) In situ images 
of the cloud just before (left) and after (right) loading into the box, 
and corresponding line-density profiles along x (bottom plots) and 
z (side plots). OD stands for optical density. The blue dashed lines 
are fits to the thermal component of the harmonically trapped gas. 
The green dashed lines are fits based on the expected profiles for a 
uniform-density gas. 
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"sheet" beams propagating along y. The green laser beams 
create a repulsive potential for the atoms and confine them to 
the cylindrical dark region depicted in red in Fig. [TJ^)- To 
create a uniform potential, we additionally cancel the gravi- 
tational force on the atoms at a 10~^ level, using a magnetic 
field gradient 1 16 |. 

As outlined in Fig.[TJb), all three trapping beams are created 
by reflecting a single Gaussian beam off a phase-imprinting 
spatial light modulator (SLM) with three superposed phase 
patterns |20|. The tube beam is an optical vortex created by 
imprinting a 24 tt phase winding on the incoming beam, the 
sheet beams are created using cylindrical-lens phase patterns, 
and the three outgoing beams are deflected in different di- 
rections using phase gradients. With a total laser power of 
Po ~ 700 mW we achieve a trap depth of ^ x 2 /iK. 

We evaporatively cool the gas in the harmonic trap down to 
T ^ 120 nK, when the cloud size is similar to the size of our 
optical box [see Fig.[TJc)] and k^T <C Vq. At this point the 
gas is partially condensed, but the BEC is lost during the trans- 
fer into the box trap, which is not perfectly adiabatic. Over 
1 s, we turn on the green light and then turn off the harmonic 
trapping, capturing > 80 % of the atoms. 

In Fig.[TJd) we show in-trap absorption images of the cloud 
just before and just after the transfer into the box trap. The 
images are taken along the y direction, using high intensity 
imaging IlllllETl with a saturation parameter ///sat ^ 150. 
For each image we show the line-density profiles along x and 

If a cylindrical box of length L and radius R is filled per- 
fectly uniformly, the density distribution along x is simply a 
top-hat function of width L. Along z, the line-of-sight inte- 
gration results in "circular" column- and line-density profiles, 
(X ^Jl — {z/Ry. In the experimental images, the edges of 
the cloud are rounded off for two reasons, both related to the 
diffraction limit of our optical setup. First, the waist of the 
532 nm trapping beams is diffraction limited to 3 /im, 
which leads to some rounding-off of the potential bottom 
near the edges of the box. Second, our imaging resolution 
is diffraction limited to ^ 5 /im, making the cloud edges 
appear more smeared out than they actually are. The green 
dashed lines in Fig. [TJd) are fits to the data based on a per- 
fectly uniform distribution convolved only with the imaging 
point- spread function. The fits describe the data well and give 
L = 63 ± 2 /im and R= 15 ± 1 /im. These values are consis- 
tent with the calculated separation of the green walls, reduced 
by the diffraction-limited wall thickness. 

After the transfer into the box trap the cloud contains N ^ 
6 X 10^ atoms at T ^ 130 nK. From this point, we cool the 
gas to below Tc by forced evaporative cooling in the box trap, 
lowering the trapping power P in an exponential ramp with a 
0.5 s time constant. Initially the trap depth is much larger than 
k^T so significant cooling occurs only for P ^ 0.5 Pq. 

Fig. |2] qualitatively illustrates the effects of evaporation and 
condensation in the box trap. We show images of the cloud 
both in situ and after t = 50 ms of TOF expansion from the 
trap. While in a harmonic trap cooling results in simultaneous 
real- space and momentum- space condensation, here it has no 




FIG. 2: Evaporation and Bose-Einstein condensation in the optical- 
box trap. Cooling is achieved by lowering the trapping laser power, 
P. We show absorption images taken after 50 ms of TOF and in situ 
[insets, with same colour scale as in Fig.[TJd)]. The bottom panels 
show cuts through the momentum distributions recorded in TOR In 
contrast to the case of a harmonic trap, no dramatic effects of cooling 
are observed in situ. However, BEC is clearly seen in the bimodality 
of the momentum distribution and the anisotropic expansion. 

dramatic effects on the in-trap atomic distribution. The den- 
sity is gradually reduced by evaporation, but the shape of the 
cloud does not reveal condensation. On the other hand, in mo- 
mentum space (i.e., in TOF) the effects of cooling are obvious 
and the signatures of condensation are qualitatively the same 
as for a harmonically trapped gas - the momentum distribution 
becomes bimodal and the BEC expands anisotropically, with 
its aspect ratio inverting in TOF. 

We now turn to a quantitative analysis of our degenerate 
quasi-uniform Bose gas. We assess the flatness of our trapping 
potential and contrast the thermodynamics of condensation in 
our system with the case of a harmonically trapped gas. 

We first analyse the momentum distribution in a non- 
condensed gas, at T ^ 110 nK. In Fig.[3ja) we plot the line 
density measured in TOF, n{z), obtained by integrating an im- 
age along X. At this point the gas is sufficiently degenerate to 
show the effects of Bose statistics. In a trapped degenerate gas 
the spatial and momentum degrees of freedom are in general 
not separable and consequently n{z) also contains informa- 
tion about the functional form of the trap. 

For simplicity let us consider an isotropic 3D trapping po- 
tential of the form V{r) oc r^; for a fixed n, potential 
anisotropy does not affect the scalings discussed below. The 
line-density distribution after long TOF then has the form 

where Qo, is the polylog function of order a = 1 + 3/n. Here 
j3 = 1/(/cbP), m is the atom mass, /i is the chemical poten- 
tial, e{z) = mz^/(2t^), and the constant 6a absorbs various 
factors such as the imaging magnification and cross section. 
Small corrections due to the initial cloud size are accounted 
for by convolving with the in-trap density distribution. 
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FIG. 3: Momentum and particle-number distribution in a quasi-uniform degenerate Bose gas. (a) TOF profile n{z) of a non-condensed cloud 
(red) and fits to the data, with a = 1 (green) and 5/2 (blue). We also show the fitting residuals, An. Inset: fitting versus a. (b) Flatness 
of the box potential: the condition S = d [log(^)] /d[log(T)] = (see text) gives a = 1.23 ± 0.03, corresponding to an effective r^^^^ 
potential, (c) The number of condensed (A^o) and thermal (N') atoms versus T. The critical temperature, Tc = 92 db 3 nK, is in agreement 
with the prediction for a uniform Bose gas. Below Tc, a power-law fit (solid red line) gives oc 7^1 '^^io.oe 



For a harmonic trap n = 2 and a = 5/2, while for a uni- 
form gas n ^ 00 and a = 1. The measured h is fitted very 
well by the gi function, with the reduced ^ 1- The ^5/2 
fit is comparatively poor: ^ 7 and the systematic patterns 
in the fitting residuals, An = h — nfit, clearly show that this 
is fundamentally a wrong functional form. Qualitatively, the 
measured momentum distribution is more "peaky" than that 
of a harmonically trapped degenerate gas. 

As shown in the inset of Fig. [SJa), the measured h can be 
fitted well using a in the range 1 — 1.7, while higher values can 
be clearly excluded. Note however that in these fits we used 
6a, /i and T as free parameters. Allowing 6a to vary gives the 
fitting function an unphysical freedom and overestimates the 
range of suitable a values. Crucially, only for the correct a 
is the best-fit value of 6a a temperature-independent constant. 
We use this fact to accurately determine the leading-order cor- 
rection to the flatness of our trapping potential. 

In Fig. |3jb) we analyse the T-dependence of the fitted 6a 
values. Here, we exclude from the fits the central 180 /im 
wide region of the cloud, which is larger than the largest ob- 
served BECs. The drift of the fitted 6a with T is described 
well by a constant slope Sa = d [log(^Q,)] /d[log(T)], and 
we see that Sa varies monotonically with a. Physically, if 
the a used for fitting is too high, as the gas cools the actual 
momentum distribution narrows faster than ga- The fit then 
increasingly underestimates T and compensates by decreas- 
ing 6a- Conversely, if a is too low the fitted 6a increases as 
the gas is cooled. Note that for each a the fitted 6a is plot- 
ted against T extracted from the same fit, and we see that the 
fitted temperatures also show the expected systematic drifts. 

From the condition Sa = we get a = 1.23 ± 0.03 and 
conclude that the leading-order correction to the flatness of 
our box potential is ex r^^^^. While we can distinguish an r^^ 
potential from a perfectly flat one, we expect this distinction 
to be irrelevant for most many -body studies. 



We now fix a = 1.23 and study the evolution of the BEC 
atom number, A^o, with T. The cylindrical BEC does not fol- 
low any simple scaling in TOF, so rather than fitting it to any 
specific shape we fit only the thermal component and simply 
count the excess atoms not accounted for by the thermal fit. 

In Fig. |3jc) we plot A^o versus T and determine Tc = 
92 ± 3 nK. From the measured atom number and box vol- 
ume, we get a consistent theoretical value for a uniform Bose 
gas, = 98 ± 10 nK. This calculation includes small finite- 
size and interaction shifts of (22H241 . about 6 % and 1 % 
respectively; the error includes the 10 % uncertainty in our ab- 
solute atom-number calibration, obtained by measuring Tc in 
a harmonic trap | 25 1. 

In Fig.jsj^c) we also plot the thermal atom number, A^^ ver- 
sus T. Below Tc, we fit the variation of A"' with a power law, 
(X T^, and get 7 = 1.73 ± 0.06. This agrees with the 
expected 7 = a + 1/2. For a saturated thermal component in 
a perfectly uniform system A^' ex T^/^, while in a harmonic 
trap N' (X T^. The thermodynamics of our gas are therefore 
very close to the textbook case of a uniform system and very 
different from the case of a harmonically trapped sample. 

Below 40 nK our temperature fits are not accurate, but we 
can cool the gas further and produce a quasi-pure BEC with 
> 10^ atoms. The 1/e lifetime of the BEC is 10 s, limited by 
the background pressure in our vacuum chamber 1 16 |. 

Finally, we verify the coherence of our uniform BEC in an 
interference experiment, using a two-trap configuration shown 
in Fig. |4ja). We add a third "sheet" laser beam to simulta- 
neously condense atoms in a small "satellite" trap, while the 
main box trap is essentially the same as in the rest of the paper. 

In Fig. |4jb) we show the interference fringes emerging as 
the two clouds overlap in TOF. Note that the interference pat- 
tern is different from the parallel fringes of constant spacing 
commonly observed in a symmetric expansion of two BECs 
1.26 J . Qualitatively, along x the main cloud expands very 
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FIG. 4: Matter-wave interference with a quasi-uniform BEC. (a) In 
situ image of a two-trap configuration, together with a sketch of the 
green-hght arrangement used to create it. (b) Interference pattern 
observed after 100 ms of TOF expansion, (c) Numerical simulation 
of the expected interference pattern. 

slowly and is still in the "near-field" regime, while the small 
BEC expands rapidly and is in the far-field regime. This leads 
to a quadratic variation of the relative phase along x and the 
reduction of the fringe period from right to left in Eig. |4jb). 
We reproduce a similar pattern in a numerical simulation of 
the Gross-Pitaevskii equation, shown in Eig.Qc). 

In conclusion, we have observed and characterised Bose- 
Einstein condensation of an atomic gas in an essentially uni- 
form potential. We expect our experiments to open up many 
new research possibilities. The box trap eliminates the need 
to rely on the local density approximation, which could be 
particularly important for studies of critical behaviour with 
diverging correlations near phase transitions. Moreover, for 
the same particle number, the quasi-uniform BEC has a sig- 
nificantly lower density (in our case ~ 2 x 10^^ cm~^) than 
a harmonically trapped one. This reduces the importance of 
three-body recombination compared to two-body interactions 
and near a Eeshbach resonance it may facilitate equilibrium 
studies of a strongly interacting, possibly unitary Bose gas 
Il27ti29ll . Dark optical traps loaded with low density BECs 
also hold promise for trapped-atom interferometry with very 
long decoherence times. Here we chose an elongated box 
shape for qualitative detection of the BEC through anisotropic 
expansion, but with our light-shaping methods traps of almost 
arbitrary geometry could be created |30|. Our methods are 
also suitable for studies of degenerate Eermi gases and are 
compatible with implementation of 3D optical lattices. 
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